Time-dependent Heavy-Quark Potential at Finite Temperature from Gauge/Gravity Duality 
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The potential between a heavy quark and an anti-quark inside the quark-gluon plasma (QGP) is studied on 
the basis of the gauge/gravity duality. A real-time complex potential Vqq (t, r) is derived from the Wilson loop 
with the Euclidean AdSs black hole metric. To make the analytic continuation from the imaginary time to the 
real time, specific variational configurations of the string world sheet in the Euclidean metric are introduced. 
Rapid approach of Vqq{t, r) to its stationary value is found at the time scale t ~ (nT)^ 1 independent of the 
't Hooft coupling A. Also, the imaginary part of Vqq (oo, r) is found to be significant above the length scale 
r = 1.72(7rT) _1 independent of A. Implications of these results to the properties of heavy quarkonia in QGP 
are briefly discussed. 
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High temperature quark-gluon plasma (QGP), whose prop- 
erties are important to understand the physics of early Uni- 
verse at 10~ 4 ~ 10~ 5 sec after the big bang, is actively 
studied in heavy-ion collision experiments at the Relativis- 
tic Heavy Ion Collider (RHIC) and the Large Hadron Collider 
(LHC). The experimental data suggest that QGP with the tem- 
perature (T) of a few hundred MeV is not a weakly interacting 
gas of quarks and gluons but rather a strongly-coupled quark- 
gluon plasma (sQGP) [1]. One of the key quantities which 
characterizes sQGP is a small ratio of shear viscosity to en- 
tropy density rj/s — 1/Att H obtained by using the holo- 
graphic duality between JV = 4 super Yang-Mills theory and 
the classical supergravity on AdSs x S 5 [3]. Another phys- 
ical quantity which can probe sQGP is the spectra of heavy 
quarkonia at finite T |4J. Such a probe was originally pro- 
posed in the context of the string breaking and Debye screen- 
ing at finite T JJ, |6[]. Later, the lattice quantum chromody- 
namics (QCD) simulations [7] and the QCD sum rule analyses 
Jit] on the basis of Bayesian technique have been attempted to 
study the melting pattern of quarkonia. It was also realized 
that the heavy-quark potential at finite T can be defined from 
the Minkowski Wilson loop in thermal medium: The potential 
thus defined is found to be time-dependent and complex by 
using thermal perturbation theory H and lattice QCD simu- 
lations iQll. (See also Ref. HO].) 



The purpose of this Letter is to derive such a time- 
dependent complex potential at finite T in sQGP on the ba- 
sis of the gauge/gravity duality in the large N c and large 't 
Hooft coupling A = gy M N c , We start with an Euclidean Wil- 
son loop We (t, r) (r (r) is the temporal (spatial) extent with 
< t < f3 — 1/T) defined in the Euclidean AdSs black hole 
metric. After evaluating We(t, r) by its gravity dual follow- 
ing the Maldacena's conjecture [12], we carry out an analytic 
continuation to obtain the Minkowski Wilson loop WM(t,f) 
with real-time t. Then, we can extract the characteristic time 
and length associated with the complex heavy-quark potential 
Vqq(£, r) obtained from Wu(t> r )- We will discuss the role 
of the imaginary part of the potential for the quarkonium melt- 



ing in our approach. We note that there is a previous attempt 
to derive a complex potential by using the gauge/gravity du- 
ality IU3I1 under Minkowski AdS^ black hole metric: We will 
mention its relation to our approach and also its limitation in 
studying the time-dependent phenomena. 

Let us start with the rectangular Wilson loop We(t, i 
the Euclidean space-time and its spectral decomposition 
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Here M is a bare heavy-quark mass taken to infinity at the 
end, and p(uj, r) is the spectral function of the heavy QQ with 
relative distance r in thermal environment. The frequency ui 
denotes energy relative to 2M, so that the spectral peaks of 
p(uj, r) associated with the bound and resonance states are lo- 
cated around u = after appropriate subtraction of the heavy- 
quark self-energy in the spatial direction. 

The time-dependent complex potential Vqq(£, r) at finite 
T in the Minkowski space-time is obtained from the analytic 
continuation of Eq.©, W M (t,r) = W E (r -> it,r), I@,[l0ll; 



V Q Q(t,r) = id t ]nW M {t,r) 
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Although p(uj,r) and We(t, r) are both real and positive, 
Wm(<, r) and the real-time potential Vqq^, r) become com- 
plex after the analytic continuation. The major questions to be 
addressed in this Letter are (i) how fast the real-time potential 
approaches to its asymptotic value as a function of t in sQGP? 
and (ii) what would be the typical time-scale and length scale 
characterizing the imaginary part of the potential? 

To evaluate We(t, r), we follow the Maldacena's conjec- 
ture 11211 and consider the gravity dual description of We (t, r) 
given by the extremum of the Nambu-Goto (NG) action in the 
background metric of the classical supergravity; 
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FIG. 1. (color online) Configuration ansatz of world sheet in the 
AdSs black hole metric with a rectangular contour on the boundary. 
Configurations I, II, and III corresponds to (a) with a > Uh, (a) or 
(b) with a — Uh, and (b) with a > uh, respectively. 




FIG. 2. (color online) Wilson loop as a function of the Euclidean 
time t for r — 1.0/nT. The solid, the dotted and the dashed lines 
correspond to the configurations I, II and III, respectively. 



where (fi = 0, . . . , 4) are the coordinates of string world 
sheet embedded in the 5-dimensional Euclidean space-time 
with a finite rectangular contour on the boundary. These co- 
ordinates satisfy the equation of motion SS-^q/SX^ = 0. We 
adopt the Euclidean AdS$ black hole metric as a dual descrip- 
tion for a strongly-coupled large iV c gauge theory at finite T; 
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where a' — I 2 with £ s being the string length and f(u) 
1 — u^/u 4 with Uh — ttTR 2 being the location of the event 
horizon in the fifth coordinate. R = (2A) 1 / 4 is a radius in the 
string unit. 

It is a formidable task to solve analytically the partial dif- 
ferential equation for string world sheet with finite rectangular 
contour on its boundary. Even if one could solve the problem 
numerically, it is not trivial to make an analytic continuation 
t — y it only with the numerical data of We (t, r) even by us- 
ing the Bayesian technique 11 1 CHI . In this Letter, we take a vari- 



ational approach to solve the string world sheet within specific 
configurations, so that one can pick up qualitative features of 
the Euclidean solution as well as its analytic continuation. The 
NG action in Eq. is proportional to the area of the world 
sheet. Therefore, taking specific configurations corresponds 
to an approximate calculation of the area. 

For small r, we consider a variational configuration I (box 
type) shown in Fig. Q2 a ) : It nas on ly one variational parame- 
ter, the bottom stringy coordinate a in the fifth direction. As r 
grows, another configuration II (wedge type) having a = Uh 
would be relevant. For r gets close to j3, we consider the vari- 
ational configuration III (anti-box type) shown in Fig.QJb). In 
the limit r — > j3, 1 (III) approaches to hanging configuration 
(two disconnected disks) of world sheet for the Polyakov loop 
correlator in Ref . Ill4ll . 

By substituting the configurations I and III into the NG ac- 
tion, we obtain a closed form with a single variational param- 
eter a. For example, the configuration I gives 
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with F being the Gauss hypergeometric function. First, sec- 
ond and third terms of Eq. © come from the area A, B and 
C in Fig. QI a X respectively. The forth term in Eq. (O corre- 
sponds to the subtraction of the quark self-energy 2uoo (t + 
r)/2n in the the temporal and spatial direction 11211 : The ac- 
tion becomes finite for —> oo after this subtraction. The 
extremum condition 5S^ G /5a = given r and r reads 




ttT 
a 



(8) 



where a scaled coordinate a = a/R 2 is introduced to elim- 
inate A. Similar equation can be derived for a(r, r) of the 
configuration III IU5I1 . 

In Fig. |2j we plot We(t, r) obtained by minimising S*ng m 
each configuration as a function of r with fixed r. We con- 
sider only the real solution of a(r, r) as physically acceptable 
one in the Euclidean space-time. It turns out that such a real 
solution of afr, r) disappear for configuration I (III) for large 
(small) t II 1611 . Although the configurations I, II and III are 
not connected smoothly, either the extended variational con- 
figurations and/or the fluctuations of string world sheet would 
eventually remove the non-smoothness. If we consider the 
branch with larger value of — S*ng given r in Fig. 12 the result 
is qualitatively consistent with that obtained from lattice QCD 
simulations 1 10] . Note that In We {t, r) behaves as ~ (^z^) 
at t — V (t — V p). However, it it is not symmetric with re- 
spect to t = (3/2 in contrast to the case of the spatial Wilson 
line correlator. 

Let us now consider the analytic continuation of We(t, r). 
Within our variational approach, it is equivalent to the analytic 
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continuation of a(r, r) by solving the master equation such as 
Eq. © with t — > it. Under the situation that the configura- 
tions I, II and III are continuous but not smooth, we primarily 
utilize the configuration I for describing the short time behav- 
ior of the potential. Furthermore, for large t and r, Vqq{t, r) 
obtained from I is found to have the same asymptotic behav- 
ior with that obtained from EH. Namely, the configuration I has 
correct information both at small t and large t. Better varia- 
tional configuration with a smooth interpolation among I, II 
and III is left for future studies. 

Under the reservation mentioned above, we solved Eq.® 
numerically with r — > it and found a unique solution a(it, r) 
consistent with a(r ~ 0, r) = \/2A/t + • • • for the configu- 
ration I. Substituting the solution into Eq. ||7), we can get the 
real-time Wilson loop Wm(£, r) as an analytically continued 
thermal Wilson loop, We(H, r). Then, by using Eq.©, the 
real-time complex potential obtained from I reads 



V Q o(t,r) 



2nr 



2ra{it, r) — {ra{it, r)) 3 + i—(ra(it, r)) 2 
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Since a(it, r) does not depend on A, the potential is propor- 
tional to y/\ irrespective of the values of T, r and t. Note here 
that a(it, r) is complex, while a(r, t) is real: This is the rea- 
son why Vnn(t, r) can become complex despite the fact that 
the Euclidean Wilson loop We(t, r) is real. 

To check if the above potential is consistent with the known 
result, we consider the asymptotic limit t — >• oo at T = in 
Eq. (0: In this case, a(it,r) — > 1/r as can be easily seen from 
Eq.©, so that we obtain, 
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2n r 
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This is the same as that given in fl2ll except that the numerical 
coefficient j- due to the box-type ansatz is smaller by 30% 



than the exact value 
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The time-dependence of the real and imaginary parts of the 
potential obtained from Eq. © are plotted in Fig.[3]for typical 
distances r = 1 / (irT) , 2 / (irT) and 3/(irT). Irrespective of 
the values of r, asymptotic values are found to be reached 
quickly near the equilibration time t cq given by 
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7rT' 



(ID 



Note that t cq is independent of A, since a and Vqq(t, r)/v2A 
are A-independent. This is in contrast to the case of weakly- 
coupled QGP (wQGP) at N c = 3 where we have t cq ~ 
10/m D ~ 10/(gT) for r = (1 ~ 3)/(ttT) (see Fig.l of 
the second reference in fl). Here the Debye mass is defined 
by m D = gTy/l + Nf /6. One finds that t cq in sQGP is not 
only parametrically different from that in wQGP, but also an 
order of magnitude shorter than that in wQGP. 

The equilibration time t oq can be regarded as a time-scale of 
the strongly coupled gluonic plasma to relax to its equilibrium 
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FIG. 3. (color online) Real and imaginary parts of the potential as a 
function of the real time t for r = 1.0/ttT, 2.0/ttT and 3.0/ttT. 



state under the external disturbance caused by a color-singlet 
heavy QQ with separation distance r at t = 0. Taking typical 
QGP temperature T = 300MeV at RHIC and LHC in Eq.lQj}, 
we have t cq ~ 0.22 fm: This is comparable to the formation 
time of the low-lying heavy quarkonia from hard-process and 
also the phenomenological formation time of QGP obtained 
from relativistic hydrodynamics Ht], so that the propagation 
of the heavy quarkonia in sQGP may well be described by 
using the complex potential in the asymptotic (static) limit. 

It is in order here to comment on the relation between our 
approach with the Euclidean AdS^ black hole metric and that 
with the Minkowski AdS^ black hole metric. Our method is 
based on the spectral decomposition of the Euclidean Wilson 
loop and its analytic continuation to the real-time on the ba- 
sis of a well-established formula in thermal field theory. On 
the other hand, if one starts from the Minkowski Wilson loop 
naively, one encounters a difficulty to select an appropriate so- 
lution 1 13]: Indeed, one needs to take into account real-time 
boundary conditions carefully [17]. One of the advantages of 
our Euclidean approach is that one can single out a unique 
configuration among various complex solutions appearing in 
the Minkowski approach. In fact, all what we need in our ap- 
proach is to consider the real stationary configuration obtained 
from the extremum of the NG action with imaginary time t, 
and then make an analytic continuation to the real-time t. 

In Figure |U we show the real and imaginary parts of the 
static potential after equilibration for different values of T. 
We normalise the potential by ttTqR 2 and plot it as a function 
of irTgr, so that the curves are A independent. Here To is an 
arbitrarily energy scale. 
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At short distance, ReVgg(oo, r) shows Coulomb-type be- 
havior given by Eq.dTOb irrespective of T. On the other hand, 
at long distance, the potential becomes deeper as T increases: 
This is consistent with what has been found in lattice QCD 
simulations where the large r behavior is dictated by twice 
of the thermal part of single-quark free energy 2Fq(T) lfl8ll . 
Note that such a deepening of the real part at long distances 
can be also seen in [13] if we use the same definition of 
ReVgg as ours by subtracting only the T-independent diver- 
gence due to the bare quark mass. We should remark here that 
the real part of our potential has a long tail ~ r -1 / 3 at large 
r and does not approach to a constant. This may be related to 
the limitation of our simple box-type ansatz. 

We find that the imaginary part of the potential shown in 
Fig- E], ImVgg(oo, r), emerges at the threshold distance 



rth 



27 



1/4 



ttT 
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Then, it grows linearly ImVqQ (oo, r ^> rth) — > 

~-^^(7rT) 2 r. This is qualitatively consistent with that ob- 
tained in fl3ll . Let us now define a "dissociation" dis- 
tance rdis where the imaginary part of the potential starts 
to dominate over the real part: For large t, Eq.® shows 
that ra is a function of ixTr only. Therefore, the condition 
ReVgg(oo,rdis) = ImVgg(oo, rais) leads to a A indepen- 
dent relation, 



'Mi: 



1.72 



(13) 



By taking typical QGP temperature T = 300 MeV at 
RHIC/LHC, we have rdi S — 0.36 fm. This is smaller than 
the Q-Q distance 0.5 fm for J/* and 0.56 fm for T(2S) but 
is larger than 0.28 fm for T(1S) lfl9ll . (For the estimate of the 
spatial sizes of heavy quarkonia, see Table 3 in |g|.) This in- 
dicates that the sequential melting of heavy quarkonia seen in 
the relativistic heavy-ion collision experiments may be closely 
related to the physics of the imaginary QQ potential rather 
than the color Debye screening. Detailed phenomenological 
studies are, however, necessary to make quantitative compar- 
ison between the theory and experiments 02OI1 . 

In this Letter, we have presented a variational derivation of 
the time-dependent complex potential between heavy quarks 
VQQ(t,r) in the strongly-coupled QGP on the basis of the 
gauge/gravity duality. Our starting point is a rectangular Wil- 
son loop in Euclidean space-time, We(t, r). Associated clas- 
sical configurations of the string world sheet under the Eu- 
clidean AdSs black hole metric were evaluated variationally 
with a single parameter a(r, r) in the fifth dimension. We de- 
rived a simple algebraic equation for a(r, r) which allows us 
to carry out the analytic continuation of We (t, r) to obtain 
the real-time Wilson loop Wu.it, f)'- Its logarithmic deriva- 
tive with respect to the real-time t is nothing but Vqq(£, r). 
Technically, the potential receives an imaginary part due to 
the process of analytic continuation. 

Resultant VQQ(t, r) is found to have interesting properties. 
It shows a characteristic equilibration time t cq ~ (7rT) -1 
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FIG. 4. (color online) The complex heavy-quark potential in the 
asymptotic limit (t —5- 00) as a function of the separation r. The 
solid and the dashed lines correspond to the real and the imaginary 
parts of the potential, respectively. To is an arbitrary energy scale. 



which is an order of magnitude smaller than the value obtained 
from the weak-coupling thermal perturbation theory. Also, the 
imaginary part of the potential starts to dominate over the real 
part at the distance r^is = 1.72(ttT) -1 which is comparable 
to the Q-Q distance of the low-lying heavy quarkonia. The 
above time scale and the length scale indicate the validity of 
using <-independent heavy quark potential in sQGP and the 
relevance of the imaginary part of the potential on the sequen- 
tial melting of heavy quarkonia in sQGP, respectively. 

We have several future directions to improve our analyses 
presented in the Letter. First of all, we need to extend our 
variational configurations so that a single smooth function for 
We(t, r) is obtained. It would be also necessary to check the 
validity of various variational ansatz by solving the equation 
of motion in Euclidean metric numerically without approxi- 
mation. Such a numerical solution may be also used to recon- 
struct the spectral function p(u>, r) by the Bayesian technique 
(e.g. the Maximum Entropy Method). 
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